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The edge states in the hybrid system of single-layer and double-layer graphene are studied in the 
tight-binding model theoretically. The edge states in one side of the interface between single-layer 
and double-layer graphene are shown to penetrate into the single-layer region when the perpendicular 
electric field is applied, while they are localized in the double-layer region without electric field. The 
edge states in another side of the interface are localized in the double-layer region independent of the 
electric field. This field-induced penetration of the edge states can be applied to switching devices. 
We also find a new type of the edge states at the boundary between single-layer and the double-layer 
graphene. 
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I. INTRODUCTION 



Recently, single-layer and double-layer graphene have 
been studied both theoretically and experimentally 1 , 
because of the interesting properties such as the 
Dirac points 2 -^, anomalous Hall effect^—, and the edge 
states 6 - - — . The double-layer graphene has attracted 
peculiar interest^ due to a band gap controlled by 
the electric field, which has been predicte d 14 ' 15 and 
observed2ii£r— . 

The edge states in the single-layer graphene and 
double-layer graphene have been studied by many au- 
thors. In the single-layer graphene the edge states ap- 
pear at the zigzag edges 6 -^ and bearded edges. If the 
system is anisotropic, the edge states also exist at the 
armchair edges 8 . The edge states in the double-layer 
graphene has been studied^ ' 10 ' 22 . The edge states at the 
interface between single-layer and double-layer graphene 
have also been studied. Transmission across the bound- 
ary has been studied theoretically using the effective- 
mass approximation 2 ^—, edge states have been stud- 
ied theoreticall y 22 ' 26 , and quantum oscillations have been 
observed in the interface 2 ^. Vacancy-induced localized 
states in the multilayer graphene has been proposed 28 . 

In this paper we study the edge states in the hybrid sys- 
tem of single-layer and double-layer graphene as shown 
in Fig. Q] We focus on the edge states localized in the 
boundary between the single-layer and the double-layer 
regions. We obtain the new edge states localized in one 
side of the interface between single-layer and double-layer 
regions with energy E ^ 0. We show an interesting prop- 
erty that the edge states at the boundary between the 
single-layer and the double-layer regions penetrate into 
the single-layer region when the electric field is applied 
perpendicular to the layers. 



II. MODEL 

We assume zigzag edges in both the first and the 
second layers. Each layer has two sublattices, which 
we call A\, B\, A2 and B2, as shown in Fig. [TJ The 
left part and the right part of the single-layer regions 
have L\l and Lm pairs of A\ and B\ sublattices in 
the x-direction. In the double-layer region there are 
L2 quartets of {A\,B\, A2,B<2) sublattices. The left 
edge in the left single-layer region has only B\ sublat- 
tice, which is labeled as n = 0, and the right edge in 
the right single-layer region (A\ sublattice) is labeled as 




FIG. 1. (color online). Single-double-single layer graphene. 
Thin black ellipses are the doublets of (Ai, B\) sites in the 
single-layer region, and thick green ellipses are the quartets of 
(Ai, B\, A2, B2) sites in the double-layer region. Thick red 
ellipses are the triplets of {A\, B\, B2) and (Ai, Bi, A2) sites 
(left and right, respectively) in the boundary between single 
layer and double layer of the a- and /3-types, respectively. 
Thick red circles are the zigzag edges. 



n = L e = Lil + L\r + L% + 3. There are two boundaries 
between single-layer and the double-layer regions. These 
two boundaries are different from each othe r 22 ' 24 " — . One 
of the boundaries has A\ , B\ and B 2 sublattices and 
the other has A\,B\ and A 2 sublattices. We call these 
boundaries as a-boundary and /3-boundary, respectively. 
The position of the a-boundary is n = L a = L\l + 1, and 
the position of the /3-boundary is n = Lp = Ln+L 2 + 2. 
Note that a-boundary and the /3-boundary always appear 
as a pair if two boundaries are parallel. We assume that 
the left boundary is the a type and the right boundary 
is the /3 type. 

We adopt the tight-binding model, where the hoppings 
between the nearest sites in the layer (A\-B\ and A 2 -B 2 ) 
are taken to be t and the interlayer hoppings between the 
nearest sites (Bi-A 2 ) are taken to be t±. We take into 
account the energy difference between layers (ei and e 2 ), 
which is controlled by the electric field perpendicular to 
the layers. 

We apply the same method which we have used in 
studying the edge states in the single-layer graphene 8 . 
Imposing the periodic boundary conditions in the y di- 
rection, we can take the wave number k y as the quantum 
number. For each k y the Hamiltonian is written as a 
L sds x L sds matrix, where L sds = 2L 1L + AL 2 + 2L 1R + 8. 
The eigenstates (vf) in the Schrodinger equation (H^f = 
E^) are vectors with L sds components (the wave func- 
tions for the first layer (^ Bl fl, ^Ai,L c , Ia^ihi an d 
^B x ,n\) where 1 < n\ < L e — 1) and the wave func- 
tions for the second layer {^B 2 ,L a , ^A 2 ,Lbi ^ 'A 2 ,n 2 > an d 
^B 2 ,n 2 , where L a + 1 < n 2 < Lp - 1)). 

In the single-layer region (1 < n < L a — 1, or Lp + 1 < 
n < L e — 1), the Schrodinger equation is written as, 

k 
_2icOS -j-^ Bl ,n-l - ttfflx.n = (E - ei)tf Al , n , (1) 

k 

-2tC08-ZV Aun+1 -m Al ,n = {E-ei)^ Bu n- (2) 

The equations in the region of double layer (L a + 1 < 
n < Lp — 1) are given by 

k 
-2tC0S^* Bli „_i -t#B u n = (E-ei)V Al ,n, (3) 

k 

-2tCOa-^^A u n+l-t^A 1 ,n 



Similarly, we obtain the equation at the right edge of 
the single-layer region to be Eq. (JTJ) with n = L e and 

-2tcos-^*B 1 ,£._i =(E-ei)*A lf L., (8) 

At the a-boundary, the equations are obtained by tak- 
ing n = L a and * A 2 x a = in Eqs. ©,((11), and @, 
since there are no A 2 sublattices at the a-boundary. 



-2tCOS^ Bl ,L a -l-t^B u L a =(E-ei)^A u L a , (9) 

-2tcofl-|*A ll i„+i-t*ii 11 i (> =(E-ei)* Bl ,L a , (10) 
~2tcos^ A2 , La+1 = (E - e 2 )V B2 , La . (11) 

The equations at the /3-boundary are obtained by tak- 
ing n = Lp and ^ B2 ,l b = 0, in Eqs. ©, (g]), and ©. 
Explicitly, the equations at the /3-boundary are given by 



k 
-2tcos-|-* Sl , L/3 _i-t* Bl , L/3 =(E-ei)-$ Al ,L 



(12) 



-2tC08-^* Al ,i 



3+1 



i* 



Ai.ifl 



(13) 

-2tCOS-^*B aii/3 _i-tj.*B 1 ,^ =(S-e 2 )*A 2 ,I,p. 

(14) 

By taking E — e±, we obtain that Eqs. (JTJ and ([2]) arc 
two independent equations for ^ Blt n and ^ Ai,n> respec- 
tively, in the single-layer regions. When ei = e 2 = 0, 
Eqs. (|3j) - (J6j) become two sets of coupled equations for 
{^A u n, ^A 2 ,n) and (^Bi,n, *s 2 ,n) in the double-layer 
region by taking E = 0. However, if ei 7^2, these equa- 
tions cannot be separated into the independent equations 
for any E. This is the origin of the field-induced pene- 
tration of the edge states into the first-layer region at the 
/3-boundary, as we will show below. 



-tx®A 2 ,n - (E - ei)* Sl ,n, (4) IIL STRICTLY LOCALIZED STATES AT k v = vr 



k 
-2*cos-|-#B 2 ,n-i -**b 2) „ 

-t^ Bl ,„ = (E-e 2 )^A 2 , n , (5) 

k 

~2t COS -^*A 2 ,n+l - t$A 2 ,n = (E- e 2 )^B 2 ,n- (6) 

At the left edge of the single-layer region we obtain 
the equation to be Eq. © with n = and ^Ai,o = 0, 
since there are no A\ sublattice at the left edge in the 
single- layer region, i.e., 

-2tcos^ Al ,i = (E-e l )^ Bu0 . (7) 



Since cosfc y /2 = at k y — it, Eqs. (JTJ - (ITU are the 
equations within the same group of n, i.e. the states at 
k y = it are strictly localized at the ellipses or circles in 
Fig. [TJ as in the single-layer graphene^. 

The energies at k y — n in the single-layer regions are 
obtained from Eqs (JTJ) and ©, as 



E.± =±t + e u 



(15) 



with (Lil + iiij)-fold degeneracy. The energies of the 
strictly localized states in the double-layer region are ob- 



tained as the eigenvalues of the matrix 



M d = 



ei -t 

-t e 1 -t x 

-t± e 2 -t 

-t e 2 



(16) 



and they are obtained to be 

ei +£2 



£, 



d.±.± 




(17) 



where Ae = e± — £2, with Z^-fold degeneracy. 

At the left and the right edges we obtain the energy as 



Er 



Er = £1. 



At the a-boundary we obtain the energies of the strictly 
localized states as the eigenvalues of the matrix 



(19) 




which are obtained as 



and 



E a ,o — € 2, 



E ay± =±t + e 1 



(20) 



(21) 



The energies of the strictly localized states at the /3- 
boundary are obtained as the eigenvalues of the matrix, 



M t 








(22) 



When |ei| <C t and |e2 1 <C t, we obtain the energies as 



Eh 



=£2 



t 2 +t\ 



tHi 



(t 2 + tl) 



2 ^4 



Ar 



(A£) 3 + 0((A£) 5 ) 



E ± =±Jt 2 +t\+ £ 2 



(2t 2 + 1\) 
2(ti+t 2 ± ) 



Ar 



± 



/ " (4 '" ?+ ^ ) .(A£) 2 + 0((Ae) 3 ). 



8(t 2 + 4) 5 / 2 



(23) 



(24) 



The eigenstates with the eigenvalues Ep^ an d Ep,± are 
obtained as 




t 1 Ae 



t 2 +*1 



0((A £ ) 2 )) 
+ 0((A £ ) 3 )) 



1 



(25) 
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FIG. 2. (color online). Energy as a function of k y for single- 
double-single graphene. There are four edge states at E — 
as shown in (b), two of them are edge states at the left and 
the right zigzag edges of single-layer regions (see Fig. [3] (a)). 
The other two edge states at E — are the edge states in the 
double-layer region at the a and the /3 boundaries (see Fig. [3] 
(b) and (c)). There exist other two edge states at E = Ep,+ 
and E = Eg t - near \k y \ = n as shown in (c) and (d). 



and 




\ 



-~^A U L 



+ 0((A£) 2 )) 
+ 0((A£) 2 )) 



(26) 



/ 



respectively. 

Since El, Er, E a $, Ep$, and Ep t ± are different from 
the energies of the macroscopically degenerate states 
(E s ,± and Ed,±,±), the eigenstates with these energies 
become the well-dehned edge states at k v ~ n, as we will 
show below. 



IV. EDGE STATES WITHOUT 
PERPENDICULAR ELECTRIC FIELD 

A. edge states with E — 

First, we study the edge states in the case of no exter- 
nal electric field (ei =£2=0). We plot the energy as a 
function of k y in Fig. [5J where we take t = 1, t± = 0.2, 
L\l = 20, L 1R = 20, L 2 = 40, and e x = e 2 = 0. As shown 
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FIG. 3. (color online). Edge states localized (a) at the left 
and the right edges, (b) at the a-boundary and (c) at the 
/^-boundary at E = and k y = 0.97T. (d) is the edge state 
localized at the /3-boundary at E (a Ep,+ or E m S/3,_ and 
k v = 0.997tt. 



in the previous section, there are four states which have 
E = when t\ = €2 = and k y = 71", i.e., El, Eh, E a ^, 
and Ep t o. Two of them are edge states localized at each 
edge in the single-layer regions (n = and n = L e ) for 
\k y \ > 27r/3, same as the single-layer graphene^. The 
edge state localized at the left edge of the single layer is 
given by 



vf 



Bi,n 



-2cos — 
2 



*; 



■1,0, 



(27) 



where < n < Ln, and other components of $ are zero. 
The edge state localized at the right edge of the single 
layer is given by 



* 



•Aljie— j 



_2cos^Y* 



A lt L e 



(28) 



where < j < Lm and other components of \& are zero. 
The other edge states with E — are localized at the 
a and /? boundaries. As shown in Appendix [XJ the edge 
state localized at the a-boundary are obtained as 



^B 1 ,L a+j 

*B 2 ,L Q +j 







-2cos ; 



* 



B 2 ,L a 



(29) 



where < j < L2 and other components of 'F are zero. 
The edge states at the j3 boundary are obtained as 



^A U L,-J 
^A 2 ,L -j 



-2cos^ 
2 



*A 2 ,L« 



1 



where < j < L2 and other components of \P are zero. 

These results are consistent with the results obtained 
in the bilayer edgei and the graphite steps 22 . The edge 
state at the /3-boundary has the finite amplitudes of the 
wave functions at A\ and A2 sites, while that at the a- 
boundary has the finite amplitude only at the B 2 sites. 

We plot the square of the absolute value of the wave 
functions in Figs. [3] (a), (b) and (c), in which we have 
taken ei = 10~ 5 and €2 = — 10~ 5 in order to lift the 
degeneracy of the edge states. We plot two edge states 
together in Fig. [3] (a), which are localized in the left and 
the right edges of the single-layer regions. There exist 
two localized states at each boundary between single- 
layer and double-layer, as shown in Fig. [3J (b) and (c). 



B. edge states with E ^ 

At the /3-boundary there exist other edge states, which 
have energy E fa Ep t ± at \k y \ « n. In Fig. |3J (d) we 
plot the amplitudes of the wave functions of the edge 
states at k y = 0.9977T and E « Ep t +. Note that if 
{^A 1 ,n 11 ^B u n 1 ,^A 2 ,n 2 -,^B 2 ,n 2 } is the eigenstate with 



energy E, {$A,,n„-* 
eigenstate with energy 
1*4 



Bi.Tii) * A 2 ,n 2 , 



-^B 2 ,n 2 } is also the 
e2 = 0. Therefore, 



P, 



|*Ai,nh 1*1?!,™, 

states with E « E : 
with E sa Ep t +. 

As seen in Fig. [5] (c 
E ~ Ep.± exist at \k y \ 



E, when e\ 

n | 2 , and |*s 2 , n | 2 for the edge 

= — -E/3.+ are the same as these 



and 



(d), the edge states at 
tt. Although the bulk ex- 
tended states with the same energy E m Ep t ± also exist 
at other values of k y , the density of states has a peak 
at that energy due to the edge states. Therefore, the 
edge states can be observed as a peak in the differential 
conductance (dl/dV) by the spatially resolving scanning 
tunneling spectroscopy (STS) 11 ' 12 at the /3-boundary. 

These edge states at E « Ep t ± can be understood 
by considering the equations in the single-layer regions 
Eqs. (HJ) and © (see Appendix [B]). 

The existence of these edge states has not been known 
before, as far as we know. Although the existence of the 
edge states at E 7^ has been suggested as a perfectly re- 
flecting states by Nakanishi et aliSi, the pure edge states 
are obtained only at E = in their paper, since they 
adopted the effective-mass scheme, which can be used 
only near the Dirac points. 



V. EDGE STATES IN THE PRESENCE OF 
PERPENDICULAR ELECTRIC FIELD 

In this section we study the edge states in the pres- 
ence of perpendicular electric field. When the electric 
field is applied perpendicular to the layers, the poten- 
tial difference between the first and the second layers, 
(ei — £2), becomes finite. Even in that case we have the 
edge states at the left and right edges in the single-layer 



(30) regions, which are given by E = ei, ^a 2 



*/ 







and ^Ai,n — (the edge states at the left edge) or 
^Bi,n = (the edge states at the right edge). The edge 
states at the a-boundary, which are given by E = e 2 and 
&Ai,n — ^Bi,n = ^ A 2 ,n = are also not affected by the 
perpendicular electric field, since the edge states at the 
a-boundary is localized only in the second layer. As seen 
in Fig.|H(a) and (b), the states with E = El = Er = e\ 
and E — E a .o — e 2 exist for \k y \ > 27r/3 even when they 
are in the upper or lower band. 

The edge states at the /3-boundary, however, is changed 
drastically by the perpendicular electric field and they 
are quite different from the edge states in the bilayer 
graphene^. 

We study the edge states at the /3-boundary with the 
energy E ss Ep t o at k y ss n in the perpendicular electric 
field in the similar method given in Appendix [Bl 

We study the single-layer region near the /3-boundary. 
We assume that the energy difference between the first 
and second layers, e± — e 2 is smaller than the hopping 
energy in the plane, t. Then we obtain 



eo| 



ei 



t 



t(e 2 - ei) 



t 2 + t 2 



<1, 



(31) 



when E « £3,0- The eigenvalues (A±) of the matrix T 
(Eq. (IB5[i l in the right single layer- region (Lp + 1 < n < 
L e — 1) are real when 



|o|<|l-|eo|| = l 



eo 



(32) 



where a = — cos(k y /2) as discussed in Appendix iBl We 
expand A± and V (defined in Eqs. (|B7|) and (|B8jl ) in a, 
and we obtain 



A. 



1 



>1, 



(33) 



(a) 



(b) 




FIG. 4. (color online). Energy of single-double-single layer 
graphene as a function of k y with different site energy in each 
layers ((a):ei = —£2 = 0.1 < t± and (b):ei = — £2 = 0.3 > 
tx-). 
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FIG. 5. (color online). Edge states at the /3-boundary in 
single-double-single layer graphene at k y — 0.9ir with ei = 
—£2 = 0.1, (a) and 0.3 (b). The localization length of the edge 
states becomes large when k y approaches to 2n/3. The edge 
states with ei = —62 = 0.1 at k y — 0.74-7T (c) and k y = 0.72n 
(d) are shown. 



Therefore, the eigenvector of T with the eigenvalue A_ is 
the edge state localized at the /3-boundary given by 






Ai*B 1|i/J+ i 



2e 
1 



Ai* Bll £ fl+1 ( l a 



(36) 



where < j < Lir. 

If Ae = 0, we obtain ^b u l p = from Eq. ((23)) an d 
we obtain ^a 1 ,l i3 +i = ^b 1 ,l^+i = from Eq. (|B24j) . 
In this case, Eq. (|36|) shows that the edge states with 
E = Epfl = and |2cos(fc y /2)| < 1 is localized only in 
the double- layer region at the /3-boundary, if Ae = 0, as 
shown Fig.[3fc). 



If Ae ^ due to the perpendicular electric field, 
^Bx,Lg and ^b ± ,Lb+i become finite for the edge states at 
the /3-boundary, resulting in the penetration of the edge 
states into the single-layer region. 

In this way the strictly localized state at the fj- 
boundary at k y = n (a = 0) and E = Ep t Q becomes 
the localized states which have the finite amplitudes 
both in the single-layer and double-layer regions when 
|2cos/e y /2| < 1 — |e |. In Fig. 0] we plot the energy as 
a function of k y in the case of the finite energy differ- 
ence between the first layer and the second layer. The 
energy of the edge states at the /3-boundary depends on 
k y as shown in Fig. 0J In Fig. [S] (a) and (b) we plot 
the wave functions of the edge states at E w Ep^ and 
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FIG. 6. (color online). Energy of single-double-single layer 
graphene as a function of k y with different site energy in each 
layers (ei = —£2 = 0.1 < t±). Two systems with different 
width {Lit = L 1R = 20, L 2 = 40 (a) and L 1L = L 1R = 40, 
1/2 = 80 (b) ) are shown. The narrower system (a) is the 
same but a close up of Fig. [4] The bulk energy gap (the gap 
between the upper blue line and lower green line at k y sa 2tt/3) 
becomes smaller as the width of the system becomes wider. 
There exists a small energy gap at k y « 27r/3 between two 
blue lines in (a), because the localization length of the edge 
state at the /3 edge becomes comparable to the width of the 
system. For the wider system (b), the energy gap between 
two blue lines becomes negligible. 



the edge states can be observed even if they are outside 
of the bulk gap. Therefore, the field-induced penetra- 
tion of the edge state into the single-layer region can be 
observed by STS. 

As shown in Fig. |4] (a), Fig. [6] (a) and (b), the edge 
states at the /3-boundary have energies in the bulk energy 
gap at \k y \ > 2tt/3 if |ei — t2 1 ^5 2£i< I n this case, the 
edge states are partially filled in the less than half-filled 
systems. Then electric current flow along the /3-boundary 
(y direction) in both single-layer region and double-layer 
region (see Fig. [S](a), (c), and (d)). The edge states are 
partially filled even when the perpendicular electric field 
is not applied. However, the contribution of the edge 
states at the /3-boundary to the electric current is only 
in the double-layer region in the absence of the perpen- 
dicular electric field, since the edge states are localized 
only in the double-layer region of the /3-boundary (see 
Fig. [3] (c)). Thus, the electric conductivity along the /3- 
boundary in the single-layer region at the /3-boundary is 
increased by the perpendicular electric field. 



VI. CONCLUSION 



k y = 0.9vt for 61 = — 62 = 0.1 and 0.3. In contract to 
the case of t\ = £2 = (Fig. [3] (c)), all components of 
the wave functions are finite in the double-layer region at 
the /3-boundary (n < Lp). The wave function has finite 
amplitudes also in the single-layer region (n > Lp + 1), 
which means the penetration of the edge state at the 
/3-boundary is induced by the electric field. When k y ap- 
proaches to 2tt/3, the energy of the edge states (the lower 
blue line) deviates from Epfi, as seen in Fig. |4] The lo- 
calization length of the edge states becomes large as k y 
approaches to 2tt/3 (see Fig. [5] (c) and (d)). The effects 
of the a-boundary (n = L a ) and the right edge (n = L e ) 
on the edge state at /3-boundary are seen in Fig. [5] (c) 
and (d). 

We mention the bulk energy gap caused by the per- 
pendicular electric field. Though the bulk energy gap 
is originated in the double-layer region, the system with 
finite width of the single-layer regions in both side as 
shown in FigfT] has the bulk energy gap, since the the 
states except for the edge states are extended in both 
single-layer and double-layer regions. When the widths 
of the single-layer regions become larger, the bulk en- 
ergy gap becomes smaller, as shown in Fig. [B] (a) and 
(b), which show the A^-dependences of the energies for 
the systems with different width. Therefore, the critical 
value of k y , at which the edge state at the /3-boundary 
has the energy between the bulk energy gap, depends on 
the width of the system. The edge states at k y = 0.74-7T 
and 0.72?T (Fig. [5] (c) and (d)) have energies between 
bulk gap (-0.0252 < E < 0.693) for the system with 
L\l = Lin = 20 and L 2 — 40, while the edge states at 
k y = 0.9tt (Fig. [5] (a)) has the energy outside of the bulk 
gap. Since STS probes the local density of state o 11 ' 12 , 



We have studied the edge states in the hybrid system 
of the single-layer and double-layer graphene. By using 
the tight-binding model, we obtain the analytic solution 
of the edge states at the boundary between single-layer 
and double-layer regions when perpendicular electric field 
is not applied. We found that the edge states at the a- 
boundary are localized only in the second layer, and the 
edge states at the /3-boundary have finite amplitude in 
both layers in the double-layer region. 

We also find the new edge states with E ss £73. ±, which 
are localized at the /3-boundary. 

When the perpendicular electric field is applied, the 
edge states at the /3-boundary are shown to change dras- 
tically. The edge states at the /3-boundary have finite 
amplitudes at all sites in both regions of the boundary. 
The penetration of the edge states induced by the elec- 
tric field can be observed experimentally by ST S 11 ' 12 . 
The bulk gap becomes smaller for the systems with the 
wider width of the single-layer regions. The edge states, 
however, are possible to be observed by STS even though 
the energy of the edge states are not located in the mid- 
dle of the bulk gap because the edge states has larger 
amplitudes of the wave functions than that for the ex- 
tended states with the same energy in the region near 
the boundary. 

We propose a simple method to observe the electric- 
field-induced penetration of the edge states. The conduc- 
tivity between two terminals placed at the single-layer 
region at the /3-boundary will become large when the 
perpendicular electric field is applied, as the edge states 
at the /3-boundary penetrate into the single-layer region. 
This electric-field-induced penetration of the edge states 
can be used as electrical devices. 



Appendix A: analytical solutions for the edge states 
with E — at the a and boundaries when ei = £2 

The equations in the double- layer region (d3j) - (JHJ) ) are 

decoupled into two groups, when ei = €2 — and E = 0. 
Eqs. Q and ([5]) are written as 

where i a + 1 < n < .Lg — 1. From this equation we 
obtain 



^ + = f-2cos^ V/ ' " W * 



ifj l)\ *B 2 ,L a J ' 



where < j < L2. Since we obtain 



(A2i 



*Bi,l„ = (-2cos^J * Bl , , (A3) 

from the equations in the single-layer region (Eq. ([1])), we 
can take $B h i Q = for |2cos(/c y /2)| < 1 and L a S> 1. 
Therefore, the edge state localized at the a-boundary are 
obtained as 







-2C0S^V *B 2 ,L 



(A4) 



where < j < L2 and other components of ^ are zero. 
In the same way we obtain from equations (|4|) and (|6|) , 



Al ' L ^ ] = |-2cos^ 



,AV^i -^iUW f 



(A5) 
where < j < L2. As in the edge state in 
the a-boundary, we can take ^Ai.l^+i — 0, when 
|2cos(A;j,/2)| < 1 and L^n ^> 1. Then we obtain from 

Eq. CE3D 



^a 1 ,l, 



tj 



t 



-* 



A2,L/3, 



(A6) 



and we obtain the edge state at the /3-boundary as 



* AuLl) - J ^ = (-2cos^) * A2 , L I < 



■j) M-'.l,.^. 



1 

where < j < L2 and other components of "J are zero. 



Appendix B: analytical solutions for the edge states 
with E ~ Ep,± at the j3 boundary 

By replacing n by n + 1 in Eq. ([1} , we can write the 
equations in the single-layer regions, Eqs. ([1} and ©, as 



e 1 \ / *Ai,«+i 
a J I * Sl ,n+l 



a 

1 e 



*Bi,n 



(Bl) 



where 



a = — 2 cos— . 
2 ' 

E-et 
e = — [ -. 



(B2) 
(B3) 



In Eq. (|B1|) . we can take 0<n<L Q orL^ + l<n< 
L e — 1, and ^a^q = ^ b x ,l q = 0. Then we obtain 



^B 1 ,n+1 J \ * Bl , 



where T is a 2 x 2 matrix given by 
/ 1 eo 



(B4) 



(B5) 



a a 



If eo 7^ 0, the matrix T is diagonalized by the matrix V 
as 



V^TV 



x+ 

A_ 



(B6) 



where A± are the eigenvalues of the matrix T, 

\ ± = ^(a 2 -el + \±^Y (B7) 



,=,^-«! + e§ + 1 _-*), ( BS, 



-2e a/D-, 



r i = l/VD-a 2 



C 



2e n 



2e 



D-a^ 



C = 2VI)(v / Z?- 



(B9) 
(BIO) 



and 



D = (a+eo + l)(a-e + l)(a+e -l)(a-e -l). (Bll) 
Note that 



A+A_ = 1. 



(B12) 



If D < 0, A_ = A+ and |A+| = |A_| = 1. In this case we 
obtain the extended states, if the boundary conditions 
at n = 0, L a , L/3, and L e are satisfied. On the other 
hand, if D > 0, two eigenvalues of T are real and either 
|A + | or |A_| is smaller than 1. In this case I^^i.l^+jI 
and \ 1 &B 1 ,L fi +j\ can decrease as |A+| J |A_p obtained from 
Eq. lB41 if the boundary conditions at n = Lp are satisfied 
by the corresponding eigenstate of T. Note that D > 
is obtained if and only if 



M|>i. 



(B13) 



Now we examine the edge states at the /3-boundary 



with E w Ep t ±. When k y w 7r and i? 
±-\/t 2 + tj_, we obtain 



£ 



|8,± 



a <1 



and 




> 1. 



(B14) 



(B15) 



Then the inequality Eq. (|B13j) is satished and the edge 
states can exist. In order to examine the edge states 
with E ~ Ep t ±, we expand D, X± and V in a when 
Id < |e | - 1, as 



D^{el-l) 2 -2(el + l)a\ 



and 



A_ 



r--,,-2 ( „!( (° ' 



po 



In this case A± are real and 

|A+|<K|A_| 



(B16) 



(B17) 



(B18) 



(B19) 



(B20) 



The eigenvector of the matrix T with the eigenvalue A+ 
is given by the first column of matrix V . Therefore, when 
^a,l,3+i and ^BXp+i (0 < j < L 1R ) satisfy the equa- 
tion, 



* 



AuLp+l 



D-a 2 + el + l 



* 



B lt L p +l 



2e 



(B21) 



we obtain from Eq. IB4I 

^AuLp+l+j \ = X jy 

«Ai* 




(B22) 



where < j < L\r. 
at the /3-boundary. 



This state is the edge state localized 



At the /3-boundary the equations for ^>A 1 ,L fl , ^B\,Lri 



*A : 

Eq. 



,is) ^AuLft + l, 

(TTJ with n = L/3 



and *Bi,l 3 

+ 1 and Eq. 



fi are 

CED, 



obtained from 



e 1 
a 






a 

1 e 




This equation is written as 



(B23) 



*Ax 



L« + l 



* 



Bi,L^ + l 



e 




/± 
















*o 


eo*_ 


a 




at 




(B24) 



This boundary condition as well as the condition that 
the wave functions should decrease exponentially in the 
double-layer region in the left part (the double-layer re- 
gion) of the /3-boundary can be satisfied by adjusting eo, 
^Ai.-La) ^B lt L0, and ^A 2 .L tt - It is indeed possible as seen 
in Fig.lSJd). These edge states at E sa -E^i are localized 
at the /3-boundary and the wave functions have the finite 
amplitudes in both sides (both single-layer and double- 
layer regions) of the /3-boundary as shown in Fig.GJJd). 



The above method can be applied to the edge states 
at the single-layer region. We consider the edge states at 
the left boundary in the left single region (n = 0) as an 
example. The boundary condition at n — is given by 



*Bi,0 




*Bi,0 



(B25) 



This state can be the eigenvector of T only when eo = 
and the eigenvalue of this eigenvector is a. Therefore, the 
edge states at the left boundary in the single-layer region 
exist only when E = e\ and \a\ — |2cos-j-| < 1, as 
obtained in the previous section. The edge states at the 
right boundary can be similarly studied by considering 
the inverse matrix of T . 
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